Differential association of key bacterial groups with diatoms and Phaeocystis spp. during spring blooms in the Southern Ocean

Abstract Interactions between phytoplankton and heterotrophic bacteria significantly influence the cycling of organic carbon in the ocean, with many of these interactions occurring at the micrometer scale. We explored potential associations between specific phytoplankton and bacteria in two size fractions, 0.8–3 µm and larger than 3 µm, at three naturally iron‐fertilized stations and one high nutrient low chlorophyll station in the Southern Ocean. The composition of phytoplankton and bacterial communities was determined by sequencing the rbcL gene and 16S rRNA gene from DNA and RNA extracts, which represent presence and potential activity, respectively. Diatoms, particularly Thalassiosira, contributed significantly to the DNA sequences in the larger size fractions, while haptophytes were dominant in the smaller size fraction. Correlation analysis between the most abundant phytoplankton and bacterial operational taxonomic units revealed strong correlations between Phaeocystis and picoeukaryotes with SAR11, SAR116, Magnetospira, and Planktomarina. In contrast, most Thalassiosira operational taxonomic units showed the highest correlations with Polaribacter, Sulfitobacteria, Erythrobacter, and Sphingobium, while Fragilariopsis, Haslea, and Thalassionema were correlated with OM60, Fluviicola, and Ulvibacter. Our in‐situ observations suggest distinct associations between phytoplankton and bacterial taxa, which could play crucial roles in nutrient cycling in the Southern Ocean.


| INTRODUCTION
The Southern Ocean (SO) plays a major role in global climate regulation through both physical and biological mechanisms, notably by the contribution of phytoplankton to the carbon cycle via the biological carbon pump and to the sulphur cycle through dimethylsulfoniopropionate (DMSP) production.Phytoplankton groups contribute differently to these processes, with diatoms being key players in the SO primary production and CO 2 drawdown (Sarthou et al., 2005) while prymnesiophytes are globally important DMSP producers (Keller et al., 1989;McParland & Levine, 2019).Besides their specific roles in elemental cycles, these groups of phytoplankton could differently affect food web structures in the SO (Hunt et al., 2021;Krumhardt et al., 2022).Phytoplankton cell size in particular influences grazing, therefore playing a key role in the transfer of cell material toward higher trophic levels and the efficiency of the biological carbon pump (Christaki et al., 2021;Eddy et al., 2021;Moline et al., 2004).For these reasons, the composition and activity of phytoplankton communities are critical determinants of SO ecosystem functioning.
The SO is the largest high nutrient low chlorophyll (HNLC) region of the world's oceans with low phytoplankton biomass during most of the year due to limitation principally by the micronutrient iron (Fe).
Spring phytoplankton blooms can, however, develop in regions with natural Fe input, such as over continental shelves or close to islands such as the Kerguelen, Crozet, and Sandwich Islands (Blain et al., 2007;Pollard et al., 2007) and are usually dominated by diatoms (Quéguiner, 2013).The prevalence of diatoms in these naturally Fe-fertilized regions is largely driven by the availability of silicic acid and diatoms' ability to outcompete other phytoplankton taxa for available iron (Marchetti et al., 2012).The onset of silica limitation and subsequent decline of the diatom blooms are usually accompanied by a phytoplankton shift from small, fast-growing to large, slowgrowing diatoms (Blain et al., 2021) to an increasing contribution of haptophytes dominated by Phaeocystis spp.(Irion et al., 2021;Salter et al., 2007).Phaeocystis is generally more abundant in coastal regions (Davidson et al., 2010) and can form extensive prolonged blooms as was reported north-west of the Ross Sea above the Australian Antarctic Ridge (AAR) (Schine et al., 2021) and north of the Crozet plateau (Poulton et al., 2007).Shifts in environmental conditions and ensuing changes in phytoplankton assemblages are expected to impact other components of the microbial community with repercussions on biogeochemical cycles and food web structure.
There is overwhelming evidence of ecological coupling between phytoplankton and heterotrophic bacteria through studies carried out both in the field and the lab.Successional changes through the course of phytoplankton blooms are mirrored in the heterotrophic bacterial community composition (see reviews by Buchan et al., 2014;Bunse & Pinhassi, 2017) that can reveal repeated patterns on multi-annual scales (Fuhrman et al., 2006;Lambert et al., 2019).These interactions may range from simple resource provision to highly specific associations encompassing a multitude of interactions from the mutualistic exchange of info-chemicals or specific nutrients to competition and antagonism (Cirri & Pohnert, 2019;Costas-Selas et al., 2024;Coyne et al., 2022).For some of these exchanges to occur, spatial proximity of the partners has long been hypothesized (Bell & Mitchell, 1972) which led to the coining of the term "phycosphere," analogous to the rhizosphere in plant-soil bacteria interactions (see (Seymour et al., 2017) for a review).In the phycosphere, high phytoplankton-derived metabolite concentrations provide a dissolved organic matter (DOM) hotspot that facilitates bacterial growth, particularly for bacteria that are chemotactic or require attachment to their eukaryotic partners.The phycosphere depends on phytoplankton cell size and exudation rates, with large (>50 µm), leaky species expected to provide a much larger phycosphere than picocyanobacteria that are estimated to fall below the theoretical lower limit of detection by chemotactic bacteria (Seymour et al., 2017).Hence, distinct phytoplankton communities with specific phycosphere properties can exert different selective pressure on coexisting bacterial taxa.
The naturally fertilized region east of Kerguelen Island has been extensively used as a natural laboratory to investigate the effect of Fe on biogeochemistry and ecosystem functioning.Pronounced and specific responses of the bacterial community during the different phases of the Fe-induced phytoplankton blooms were observed (Hernandez-Magana et al., 2021;Landa et al., 2016;Obernosterer et al., 2011;West et al., 2008).However, these previous studies focused mostly on free-living bacteria (0.2-0.8 µm).The major aim of the present study was to identify associations between phytoplankton and particle-attached bacteria by focusing on two size fractions (0.8-3 and >3 µm) that should separate the larger and smaller size classes of phytoplankton (Biggs et al., 2019) and the results presented in this study complement observations in the freeliving fraction from the same period and study sites (Dinasquet et al., 2022).To achieve this objective, heterotrophic bacteria and phytoplankton diversity and activity were assessed in spring phytoplankton blooms by metabarcoding of the 16S rRNA gene for bacteria and the RubisCO rbcL gene for phytoplankton.As a proxy for activity (protein synthesizing potential; (Blazewicz et al., 2013), cDNA fragments generated from RNA from each marker gene were also analyzed by the same approach.

| Study area
The samples for the present study were collected in the naturally Fe-fertilized and high nutrient low chlorophyll (HNLC) waters off Kerguelen Island during the KEOPS2 (Kerguelen Ocean and Plateau Study 2) cruise aboard the French R/V Marion Dufresne II ship in the austral spring (10 October to 20 November 2011) (Figure 1).

| Sample collection
Seawater samples were collected from the four stations at 3-4 discrete depths using a CTD equipped with 12L Niskin bottles (General Oceanics) (See Table 1).One set of samples from the upper 3 depths (20-150/160 m) were size fractionated by prefiltering 9L of seawater through a 65 µm mesh and then filtering onto 10 µm poresize polycarbonate filters (PC, Nuclepore, 47 mm diameter).
The second set of samples was collected from 4 depths (20-300 m), and passed through a 25 µm mesh, before sequentially filtering onto 3 µm and 0.8 µm PC filters followed by a 0.2 µm Sterivex cartridge filter.The volumes filtered ranged between 5 and 7 L except for the 20 and 300 m depths at station R where 2.5 and 3 L were filtered respectively.All filters were stored frozen (−80°C) until further analysis in the laboratory.In this study, the samples from the 10 , 3, and 0.8 µm filters were analyzed whereas the 0.2 µm Sterivex cartridge filters were analyzed separately (Landa et al., 2016).

| DNA and RNA extraction
DNA and RNA were extracted simultaneously from the filters as described previously (West et al., 2016).Briefly, bacterial lysis was achieved by adding to the filters 425 µL lysis buffer (40 mM EDTA, 50 mM Tris, 0.75 M sucrose) and carrying out three freeze-thaw cycles (liquid nitrogen −65°C) followed by lysozyme treatment (final concentration 1 mg/mL) at 37°C for 45 min.Proteinase K and SDS were added (final concentration 0.2 mg/mL and 1% respectively) to the filters and incubated at 55°C for 1 h.Purification of DNA and RNA was carried out with the Qiagen AllPrep DNA/RNA extraction kit using 1.55 mL RLT + buffer containing β-mercaptoethanol according to the manufacturer's instructions.The quality of the DNA and RNA was verified by agarose gel electrophoresis and quantified by Picogreen and Ribogreen (Invitrogen) respectively.

| cDNA synthesis
cDNA was synthesized from RNA immediately after extraction using the Superscript VILO kit (Invitrogen) according to the manufacturer's instructions using 100 ng RNA.

| Illumina sequencing
To determine bacterial diversity and activity, bacterial 16 S rRNA gene fragments of about 450 bp were amplified from 36 DNA samples and 35 cDNA samples (see Table A1) using the primer pair 341F (CCTACGGGNGGCWGCAG) and 805R (GACTACHVGGGTATC-TAATCC) described previously (Klindworth et al., 2013).Phytoplankton diversity and activity were assessed by amplifying rbcL fragments of 554 bp from 16 surface DNA samples and 4 surface cDNA samples from the 2 bloom stations A3 and F-L (Table A1) using the Form ID RbcL gene primer set (F: GATGATGARAAYATTAACTC, R: ATTTGDCCACAGTGDATACCA) described previously (Wawrik et al., 2003).The 341F and the rbcL forward primers were tagged at the 5ʹ end with different 7 bp tags for each sample, and that were chosen from a set of tags designed to be robust to substitution, deletion and insertion errors incurred in massively parallel sequencing (Faircloth & Glenn, 2012).We also included in parallel a control con-  A1) to give a combined fraction size of >3 µm <65 µm (pooled at equal concentrations of DNA/cDNA).
We cannot rule out that phytoplankton falling between 10 and 25 µm may be overestimated due to the overlap between these 2 size fractions.The samples from the 0. pipeline MiSeq Control Software 2.4.1.3,RTA 1.18.54.0,The error rate was measured by spiking the library with about 0.5% of a PhiX library and mapping the reads onto the PhiX reference genome.

| Sequence preprocessing
Paired-end Illumina sequence reads were initially processed using version 8 or version 9 of the USEARCH-64 package (Drive5) for read merging, quality filtering, parsing of the reads according to the project and primer removal, with some commands from QIIME version 1.9.1 (Caporaso et al., 2010) and mothur version 1.33.3 (Schloss et al., 2009).
Several criteria were chosen in the preprocessing steps to minimize errors and reduce over-inflation of diversity (Bokulich et al., 2013); (a) only three mismatches were allowed in the overlapping region when merging the paired-end reads, (b) quality filtering was carried out after merging using a stringent expected error of 1.0 (Edgar & Flyvbjerg, 2015), (c) zero mismatches were allowed in the barcode when demultiplexing, and (d) exact matches to both primers required.More details are found in Appendix 3.

| Definition of zOTUs
The UNOISE3 algorithm from the USEARCH-32 package (version 11) was used to denoise the 16S rRNA or rbcL reads (removal of sequencing errors and chimaeras), allowing the recovery of correct biological sequences, also termed zero-radius OTUs (zOTUs) (see Edgar, 2016) for the description and validation of the previous algorithm UNOISE2).

| Classification of 16S rRNA sequences
zOTUs were classified with the USEARCH-32 v11 classifier SINTAX and the SILVA v123 sequence database using a cutoff of 0.8 (see Appendix 3).zOTUs were aligned using the mothur SILVA SEED T A B L E 1 Brief description of the study sites.For bacterial diversity, the PhyloSeq object was filtered to remove Archaea, Eukaryotic, mitochondrial, and plastid sequences.
Of the 232 zOTUs identified as plastids, around 93% were only classified as far as Order level as Chloroplasts using the SILVA 123 database.To obtain a classification at a lower taxonomic level, all zOTUs classified as Chloroplasts were extracted and classified further with the SINTAX classifier as carried out for the bacterial sequences, using the Phytoref Database (Decelle et al., 2015;see Appendix 3) before reimporting into Phyloseq.

| Taxonomic assignment of rbcL sequences
In the absence of an rbcL sequence database for marine phytoplankton, zOTU sequences were analyzed by NCBI Blastn (Altschul, 1997) and the hit with the highest percentage similarity including a named species was used to assign the taxonomy with the taxonomic levels retrieved from the WoRMS database (WoRMS Editorial Board, 2016).Sequences with >99% sequence identity to a named species were assigned to species level unless two different species showed the same % identity.In this case, the zOTU was assigned as far as the genus level.zOTUs showing a % ID of 97% or greater were assigned to the genus level.Blastn results with hits <97% were assigned to the lowest taxonomic level possible coherent with the different hits proposed in the results table (e.g., for Zotu74, there were only two hits at 96% with the named species Teleaulax sp. and Phagioselmis sp, both belonging to the family Geminigeraceae and therefore the zOTU was assigned at this level).In the case of blastn results that retrieved only uncultured sequence hits, even after allowing 250 hits, uncultured clones that had been situated in phylogenetic trees in published data (e.g., Bhadury & Ward, 2009) were used to attempt to assign the sequences at least at the class level.
zOTUS assigned to non-marine groups including the Phylum Streptophyta (land plants) and the Class Trebouxiphyceae as defined by the World Register of Marine Species (WoRMS) were filtered out of the OTU table as most likely representing cross-talk OTUs from the lichen project analyzed in the same run.
The rbcL sequences were imported into ARB (arb-6.0),aligned and a neighbor-joining tree was constructed.The tree was exported in Newick format to be combined with the zOTU table, tax_table, and meta-data using the Phyloseq package.

| Data analysis
Data manipulation and analysis was performed using various packages of the R (version 3.5.1)platform (ape, phyloseq, vegan, mixOmics, ggplot2, dplyr, tidyr, factoextra, dendextend, matrixStats, stringr, reshape2, tidyverse).Alpha diversity metrics were calculated with the phyloseq package and the Shannon Entropy values (x) were transformed to effective species diversities according to (Jost, 2006) by taking the exponential of the values (exp(x)).Rarefaction of zOTU tables was chosen before clustering methods since this normalization method may better reflect the clustering of samples according to their biological origin (Weiss et al., 2017).Nonmetric dimensional scaling (NMDS) plots were generated by Phyloseq using Bray-Curtis dissimilarity (upper three depths) and dendrograms were constructed using hclust() and Ward.D2 clustering.The significance of the bacterial community structure differences represented in the nMDS plots was tested by PERMANOVA using the adonis2 routine in vegan after first verifying the homogeneity of group dispersions with betadisper.
For the figures focusing on specific taxonomic groups (Figures 4-6, A6 and A7), a reduced data set was used to eliminate potentially spurious bacterial zOTUs from Illumina cross-talk artefacts.In brief, the variance was calculated for the zOTUs as carried out previously (West et al., 2018), retaining only those zOTUs with a variance >1 × 10 −7 .This reduced the zOTU table to 679 taxa but still retained 91.3% ± 2.3% of the total read count when considering the upper three depths of the four stations or 89.7% ± 4.3% when including the 300 m depth samples.(Figure 1).All stations were located south of the polar front except for station F, and their hydrographic conditions are described in detail in (Park et al., 2014) and(d'Ovidio et al., 2015).Station R was sampled as a reference HNLC site west of the Kerguelen plateau with surface waters characterized by low concentrations of dissolved iron (dFe; <0.088 nM) (Quéroué et al., 2015) and Chl a (<0.32 µg L −1 ) (Lasbleiz et al., 2016).The three stations subject to natural Fe fertilization had overall higher concentrations of dFe (0.18-0.66 nM) in the wind mixed layer (Z ML ) and Chl a concentration varied considerably, with values ranging from 0.9 (E5) to 4.0 µg L -1 (F) among the Fe-fertilized stations.This patchiness in Chl a and other biological parameters was due to spatial and temporal variability in the bloom development (Lasbleiz et al., 2016).By contrast, the major inorganic nutrients N and P were similar across sites and characteristic of this region (Blain et al., 2015).Stations R and A3 had deeper Z ML (105 and 153 m, respectively) than stations F and E5 (38 and 46 m, respectively) (Table 1) with the euphotic layer depth (1% surface PAR) roughly equaling the Z ML at all sites except for station A3, where PAR penetrated to only a small fraction of the Z ML .species in the Southern Ocean (e.g., Wright et al., 2010), but would not target the green algal groups belonging to Chlorophyta and cyanobacteria or dinoflagellates (e.g., Tabita et al., 2008).We chose to use rbcL since it is plastid-encoded and therefore excludes heterotrophic protists that are often overrepresented in 18S rRNA sequence libraries (Georges et al., 2014;Heywood et al., 2011).The rbcL gene also performs better than the 18S rRNA gene to distinguish closely related taxa (Evans et al., 2007).In addition, since rbcL is a key gene involved in photosynthetic carbon fixation, the analysis of rbcL RNA transcripts could be used as a proxy for phytoplankton activity (Wawrik et al., 2003).
After paired-end read joining and quality filtering, 0.9 million 16S rRNA gene reads were denoised followed by clustering into 1719 bacterial zOTUs (100% identity; zero-radius OTUs) and 232 plastid zOTUs (see Experimental Procedures).Similarly, just over 0.5 million rbcL reads were clustered into 222 rbcL zOTUs.Although 16S rRNA zOTUs assigned to plastids are often filtered out in bacterial diversity studies, here they were retained as a useful comparison to the picture of phytoplankton diversity revealed by the rbcL sequences.The proportional abundance of the 16S rRNA bacterial sequences was highest in the 0.8-3 µm fraction whereas plastid sequences were proportionally more abundant in the >3 µm size fraction except for the 300 m samples where bacterial sequences dominated (Figure A1).To examine the diversity and composition of the active fraction of bacterial and phytoplankton communities, 16S rRNA and rbcL gene sequences were also generated from reverse-transcribed RNA (cDNA) for the majority of the samples for 16S rRNA and from four surface samples at the A3 and F bloom samples for rbcL (see Table A1 for details).
The alpha diversity index Shannon Entropy was calculated from rarefied zOTU tables and then transformed to true diversities or "effective species number" for both phytoplankton and bacteria (see Experimental Procedures).This is defined as the number of evenly represented species to give the value of the calculated index which allows an intuitive comparison between the samples (Jost, 2006).For phytoplankton, the 16 S rRNA and rbcL data gave diversity estimates of a similar magnitude for a given sample (Table A2).The R and E5 >3 µm size fraction phytoplankton diversity was over threefold higher than that at the bloom stations (A3 and F) whereas the 0.8-3 µm fractions of all stations gave similar true diversities of around 14-20.
Similarly, microscopic observations revealed a higher diatom diversity, corresponding mainly to the >3 µm size fraction, at the E5 and R sites as compared to A3 and F (Lasbleiz et al., 2016).Bacterial diversity was generally higher in the >3 µm fraction than in the 0.8-3 µm fraction in the two surface depths.When ranks of bacterial diversity and depth were compared (1-4, lowest to highest diversity or surface to deep; Table 2), diversity was strongly positively correlated with depth (R 2 = 1) for the 0.8-3 µm fraction, while for the >3 µm fraction there was no significant correlation between diversity and depth.WEST ET AL.
| 7 of 32     In the >3 µm fraction of the upper two surface depths, Bacillariophyceae accounted for >90% of the sequences at stations A3 and F but exhibited lower relative abundances (>60%) in the same size fraction at the E5 and R stations (Figure A3).The lower relative abundance of diatoms at E5 and R compared to A3 and F was also observed by (Lasbleiz et al., 2016) who assessed the abundance of different phytoplankton groups by microscopy counts and flow cytometry.Diatoms accounted for 72%-74% of the phytoplankton groups in the euphotic zone at A3 and F but only 38% and 7% respectively at E5 and R (Lasbleiz et al., 2016) A1), the two analyses yielded consistent relative abundances of the two major classes Bacillariophyceae (diatoms) and Prymnesiophyceae (haptophytes).Notable differences between the two marker genes included the higher relative abundance of Pelagophyceae and Cryptophyceae zOTUs in the 16S rRNA data set compared to the rbcL data set, whereas the silica-containing Bolidophyceae group that is closely related to diatoms was observed in the rbcL but not the 16S rRNA data set.It is not uncommon to see differences in the relative abundance of phytoplankton groups when using different primer pairs due to the effects of primer bias as observed by (Shi et al., 2011).In their study, Pelagophyceae and Dictophyceae sequences were favored with one pair whereas Prymnesiophyceae dominated the sequences obtained with the second pair, underlining the utility of multiple marker gene comparisons to obtain a fuller picture of phytoplankton diversity.
At the Order level, in the >3 µm fraction, both marker genes indicated a dominance of Thalassiosirales diatoms in the bloom stations A3 (73% of total sequences) and particularly at F, where these sequences accounted for >90% of the total rbcL sequences (Figure 2).Surprisingly, at A3, higher relative abundances of the diatom order Chaetocerotales were observed in the rbcL sequences of the smaller 0.8-3 µm fraction compared to the >3 µm fraction whereas the converse was true at R. For the 16S rRNA sequences, the A3 and F stations shared the same main zOTU that was assigned at Order level (Thalassiosirales; data not shown) whereas the rbcL sequences distinguished the dominant zOTU at these two stations (Figure 3).A single rbcL zOTU assigned as Thalassiosira ritscheri accounted for 60%-75% of the sequences at F and less than 0.4% at A3, whereas the second most dominant zOTU was assigned as Thalassiosira sp.
(97% similarity to Thalassiosira antarctica) and accounted for 49%-54% of the sequences at A3 and less than 2.1% at F. The >3 µm phytoplankton communities at the HNLC station R and E5 showed a more diverse and even distribution with higher relative abundances of Fragilariopsis spp, Pseudo-nitzschia spp., Chaetoceros spp.and Eucampia spp.The haptophyte rbcL sequences recovered from the 0.8-3 µm fraction were dominated by two roughly equally abundant Phaeocystis zOTUs at the bloom stations A3 and F (zOTU3 and zOTU6) and a single Phaeocystis antarctica (zOTU3) at E5 and R. The major Thalassiosira zOTU1 also accounted for a significant fraction of the rbcL sequences in the 0.8-3 µm fraction at F (Figure 3).
We compared the relative abundance of diatoms as assessed by rbcL gene sequencing from the >3 µm samples to microscopy counts carried out on CTD rosette samples (Lasbleiz et al., 2016 and Table A3).To allow comparisons, the species identified by microscopy and the rbcL zOTUs were grouped, when possible, at the order level, apart from a few low-abundance species that did not fall into the orders listed and were grouped instead under "Class_Bacillariophyceae."The main differences were observed for the Thalassiosirales and Chaetocerotales at the bloom stations where our rbcL sequence data set appeared to overestimate Thalassiosirales abundance: in particular the dominance of Thalassiosira spp. at A3 diverges from the microscopy counts where Chaetoceros (Hyalochaete) spp.
accounted for 60% and 80% of carbon biomass and live cell abundance respectively (Lasbleiz et al., 2016, Table A3).This could be explained by our 65 µm prefiltering step that could have reduced the number of the Chaetoceros chain-forming diatoms in the sequence datasets or less efficient lysis of certain species with highly silicified frustules at the nucleic extraction step (Luddington et al., 2016).For station R, the diatom order abundances were relatively similar between microscopy counts and sequencing data except for the absence of Thalassionematales in the rbcL data.
The analysis of the comparatively small RNA (cDNA) rbcL data set obtained from the surface depth of 20 m at the bloom stations A3 and F provided further evidence that the dominant diatoms affiliated to Thalassiosira spp.were also active, given the similar relative abundances in the >3 µm fraction for RNA or DNA.Interestingly the 0.8-3 µm size fraction showed more pronounced differences between the RNA and DNA data sets, notably at station A3.Whereas Thalassiosira sequences were also the major contributor to the station F RNA sequences, at A3, Chaetoceros spp.accounted for >65% of the RNA sequences (Figure A4).
Therefore, despite the apparent overestimation of Thalassiosira spp. in the rbcL DNA sequence data, we were nonetheless able to detect the main diatom species and highlight a dominant Chaetocerales activity in the RNA data set at A3 in agreement with (Lasbleiz et al., 2016).Furthermore, we were also able to assess the presence, respectively, have been applied to marine and freshwater bacterial communities (Denef et al., 2016;Hunt et al., 2013).16S rRNA gene sequences from RNA and DNA were also compared with metatranscriptomic analysis on the samples to investigate the impact of a phytoplankton bloom on bacterial diversity (Wemheuer et al., 2014).
The high coupling between the composition of the total (DNA) and the active (RNA) bacterial communities observed only in bloom stations suggests a stronger role of bottom-up control in the shaping of community composition in response to an increased supply of photosynthate.In contrast, higher fractions of dormant cells and slow-growing bacteria could further explain the distinct clustering between the total and active bacterial communities at the HNLC site.
Bacterial community composition was significantly different between sample stations, size fractions, sample type and depth, as illustrated in the nMDS plots in Figure 4 and supported by the PERMANOVA analysis in Table 2.These differences can also be observed by the relative abundance of the different taxonomic groups in Figure A6 and the hierarchical clustering shown in Figure A5.Communities in both size fractions from below the Z ML clustered separately from those within the Z ML and this was particularly noticeable at station F where the Z ML was relatively shallow (See Z ML in Table 1 and Figure A5).At a glance, even at the order level, there were clear differences in the taxonomic composition of the different size fractions.Our results are in line with studies from other marine environments showing that particles of different sizes harbor distinct bacterial communities, and that specific bacterial taxa are found preferentially associated with certain particle sizes Mirror column chart to compare the relative mean abundance (upper two surface depths) of the 40 most abundant rbcL zOTUs in the 0.8-3 µm and >3 µm fractions at the four stations.(Mestre et al., 2017(Mestre et al., , 2020;;Yung et al., 2016).To reveal differences in the preferences of different bacterial taxa for the particulate niches represented by the two size fractions, a Particle-Association Niche Index (PAN-Index) was calculated for each zOTU as described previously (Mestre et al., 2018;Salazar et al., 2015).The PAN-Index of taxa found predominantly in the 0.8-3 µm fraction showed a value towards 0 whereas taxa found more often in the >3 µm fraction had a value closer to 1 (Figure 5).
In the mixed layer depths (Z ML ), the >3 µm particle niche was Rhodobacteraceae.Synechococcus zOTUs were also most abundant in the 0.8-3 µm fraction as observed previously (Mestre et al., 2017), most likely due to their larger cell size than heterotrophic bacteria.
However, relatively high abundances (5.5%-6%) were also observed in the >3 µm fractions at stations R and F, respectively, which could be explained by the formation of aggregates, attachment to particles or due to symbiotic relationships.Synechococcus-like symbionts have been detected in dinoflagellates, tintinnids and radiolarians (Foster et al., 2006) and in Foraminifera (Bird et al., 2017).Furthermore, high abundances of Synechococcus psbO marker gene sequences were also recovered in metagenomic libraries from large-size fractions collected by Tara Oceans (Pierella Karlusich et al., 2023).

| Abiotic versus biotic variables influencing bacterial community structure
To explore the influence of abiotic versus biotic variables on the bacterial community structure, partial Mantel tests were performed between dissimilarity matrices of bacterial 16S rRNA gene zOTUs versus phytoplankton zOTUs assessed either by plastid 16S rRNA gene or by rbcL gene diversity (biotic) and versus an environmental variable matrix (abiotic).In the upper two surface depths (20 m and 60-80 m) bacterial communities were significantly correlated with phytoplankton (based on 16S rRNA) for both size fractions, but not to environmental variables (Table A4).When all depth layers were considered, bacterial communities were correlated with phytoplankton but more strongly with environmental variables, highlighting the impact of the environmental conditions set by depth on bacterial community composition.For the rbcL data, only the 20 m and 60-80 m surface samples were available, and for both size fractions, bacterial communities showed significant correlations with phytoplankton but not with environmental variables.These observations The similarity of the bacterial communities in different size fractions collected from the four stations in three depth layers (20, 60-80, and 150 m) assessed by 16S rRNA marker gene analysis of DNA or RNA (cDNA) using Bray-Curtis distances and nMDS.The same nMDS is shown in two separate panels to show the separation of the samples by fraction and station (left) and by depth (right).Size fraction 0.8_ micron (0.8-3 µm), 3_micron (>3 µm).
are in line with the conclusions from several recent studies showing that biotic interactions could be more important for structuring plankton community composition than abiotic factors (Arandia-Gorostidi et al., 2022;Lima-Mendez et al., 2015;Liu et al., 2019;Needham et al., 2018).

| Phytoplankton-bacteria co-occurrences
Given the strong correlation between the phytoplankton and bacterial communities indicated by the partial Mantel tests, we identified specific correlations between the most relatively abundant phytoplankton and bacterial zOTUs (representing >1% of the total sequences in at least one sample) for the upper two surface depths (20 m and 60-80 m) by sparse partial least squares (PLS) regression provided by the MixOmics R package.Correlations were done between the rbcL zOTUs and 16S rRNA bacterial zOTUs from DNA samples (Figure 6a) and also between the 16S rRNA plastid and bacterial zOTUs from DNA and from RNA (cDNA) potentially representing the active zOTU interactions (Figure A7).Further information on the relative abundance of the strongly correlating bacterial zOTUs and their potential activity is presented below the heatmap (6b).
This correlation analysis presented in Figure 6a identified 56 and 43 strongly positively correlated (r > 0.6) phytoplankton and bacterial The size fraction preference of zOTUs in the mixed layer depths (Z ML ) at the four stations inferred by calculating a Particle associated Niche Index (PAN index) and illustrated in a boxplot with zOTUs grouped at the family level.zOTUs with a PAN index close to zero are prevalent in the 0.8-3 µm fraction whereas those with a PAN index toward one are more abundant in the >3 µm fraction.Subsection I, Family I corresponds to Synechococcus sp., Family Incertae Sedis to Marinicella sp. and the family Dietziaceae is affiliated to the Corynebacteriales order.Families are color-coded according to order as illustrated in the legend.

(a)
F I G U R E 6 (a) Clustered image heatmap showing the most significant correlations between the surface bacterial 16S rRNA zOTUs and phytoplankton rbcL zOTUs (upper 2 surface depths).The main bacterial zOTU clusters are outlined in blue and illustrated similarly in (b).(b) For clarity, the summed relative abundances (DNA) of bacterial surface zOTUs in the 2 size fractions (0.8-3 and >3 µm fraction) and their relative activity (ratio of RNA:DNA relative abundance) is shown only for the two main bloom stations A3 and F (omitting the meander station E5) for comparison with the HNLC station R.

(b)
F I G U R E 6 (Continued).
WEST ET AL.
| 13 of 32     zOTUs, respectively.Together, the phytoplankton zOTUs correlating to bacterial zOTUs contributed from 59.8% to 89.9% (n = 12) of total phytoplankton communities at the two surface depths when excluding station F, for which the relative abundances of correlating zOTUs was significantly lower (8.9% to 63.7%, n = 4), in particular in the >3 µm fraction (Table A5).In contrast, the contribution of bacterial zOTUs correlating to phytoplankton was lower but more stable (31.0%-49.9%n = 16) with no particular patterns between stations (Table A5).The difference between phytoplankton and bacteria in the contribution of correlated zOTUs to total communities could reflect a general tendency of phytoplankton to associate with specific bacteria while from the bacterial point of view, a significant proportion of the community could behave like generalists and thrive on particles, irrespective of the identity of phytoplankton present within the particles.This could suggest that the bacterial zOTUs identified as correlated to phytoplankton zOTUs could be involved in true partnerships with specific phytoplankton taxa.
Bacterial and phytoplankton zOTUs with strong positive correlations (r > 0.6) were distributed in several distinct clusters driven in part by the compositional differences observed between the two size fractions (Figure 6a).Cluster 1 contained taxa associated with the 0.8-3 µm fraction with the strongest correlations between the Prymnesiophyceae and Pelagophyceae phytoplankton zOTUS with SAR11, SAR116, and OCS116 clades, Magnetospira and 3 Planktomarina (Roseobacter RCA clade) zOTUS.These latter zOTUs (zOTU48, 93, 34) were the only members of the cluster showing a higher relative activity as observed in Figure 6b (excluding Synechococcus).In this cluster, the presence of typical pelagic oligotrophs such as SAR11 and the Roseobacter RCA clade (Planktomarina) and their high correlations with Pelagophyte (Pelagomonas sp.) and Haptophyte (Phaeocystis sp.) zOTUs suggest that for some of these bacterial taxa, interactions with phytoplankton might be more prevalent than previously recognized.Such interactions could rely on the exchange of organic sulfur compounds such as DMSP or sulfonates, which prymnesiophytes are known to produce, that could fulfill the organic sulfur requirements of taxa that have lost the sulphate assimilation pathway such as the SAR11 (Durham et al., 2019;Tripp et al., 2008).Our data suggest that physical association, either transient or more stable, could play a part in the nutrient exchange between these organisms, as was also suggested by a few other studies that identified SAR11 ecotypes in >0.8 µm size fractions (Allen et al., 2012;Mestre et al., 2017).Our findings differ from those of a Phaeocystis-dominated bloom in an Antarctic polyna, where the Phaeocystis cells were in a colony form, and consequently fell into the larger >3 µm fraction, correlating instead with the Gammaproteobacteria SAR92 (Delmont et al., 2014).In our study, SAR92 was most abundant at A3 but did not exceed more than 2.8% of the total bacterial community and did not show a particular size fraction preference (Figure 5b, Porticoccaceae), although the correlation with one SAR92 zOTU was strongest with diatom zOTUs (Figure 6, cluster 5).In a previous study, samples were taken during the peak and decline of the diatom bloom at A3, and SAR92 was both relatively abundant and active in the surface waters (Obernosterer et al., 2011;West et al., 2008), supporting the correlation observed in this study.
Cluster 5 highlighted correlations between several Thalassiosira, Pseudo-nitzschia, Bacterosira, and Chaetoceros zOTUs, with two Roseobacter zOTUs (Sulfitobacteria and Erythrobacter), Bacteroidetes zOTUs (Polaribacter, Flavobacterium and Owenweeksia), SAR92, and Sphingobium.These zOTUs were generally more abundant at A3 than F or R and also showed higher relative activity, particularly for the Roseobacters and Sphingobium (Figure 6b).Sphingobium is often described as a polycyclic aromatic hydrocarbon degrader but reports of the interactions of this taxon with phytoplankton are scarce and would warrant further investigation.Roseobacter OTUs were also amongst the most active members of the bacterial communities dominating the peak of the bloom at the same station A3 sampled during the KEOPS1 study (Obernosterer et al., 2011;West et al., 2008).Associations between diatoms and members of the highly diverse Roseobacter clade are among the best-characterised marine microbial partnerships that have important repercussions for marine biogeochemical cycles (e.g., Buchan et al., 2014).Many phytoplankton are vitamin B 12 auxotrophs (Croft et al., 2005) and therefore these associations typically revolve around the trading of essential metabolic currencies, such as the provision of the sulfonates in exchange for vitamin B 12 , known to be synthesized by many Roseobacter clade members (Durham et al., 2015).There are several reports of the specific association of Sulfitobacter strains with Pseudonitzschia species that promote the growth of the diatoms through their provision of B vitamins (Andrew et al., 2022) or by the secretion of the growth-promoting hormone indole-3-acetic acid (Amin et al., 2015).An extensive study of Roseobacter biogeography from an Australian time series revealed the co-occurrence of Planktomarina and Sulfitobacter with high DMSP-producing phytoplankton, such as Phaeocystis (O'Brien et al., 2022).Although also correlated to Thalassiosira, one Sulfitobacter amplicon sequence variant (ASV) showed the highest correlations to Phaeocystis that exceeded those of SAR11 and SAR86 (O'Brien et al., 2022).This is in contrast to our size fractionated data where the Planktomarina and Sulfitobacteria were significantly correlated to Phaeocystis and diatoms respectively.For bacterial zOTUs preferentially found in the >3 µm size fraction and correlating to diatom zOTUs, the majority belonged to the phylogenetic groups Flavobacteria, Gammaproteobacteria, and Roseobacter that are typically found in association with phytoplankton (Abell & Bowman, 2005;Liu et al., 2019;Tran et al., 2023).The distinct associations observed for specific diatoms could result from substrate preferences that vary among bacterial taxa (Bunse et al., 2021;Krüger et al., 2019;Xing et al., 2015) responding to differences in the composition of organic matter released by diatom species (Mühlenbruch et al., 2018).Genomic analyses of Flavobacteria demonstrate that they possess many genes involved in adhesion and gliding motility as well as a unique secretion system (Fernández-Gómez et al., 2013;Gavriilidou et al., 2020) that allows them to employ their diverse repertoire of enzymatic capabilities for the degradation of complex organic matter (Kappelmann et al., 2019).These traits are thought to facilitate a particle-attached lifestyle, including physical association with phytoplankton.Our observation of several Flavobacteria members with relative abundances correlating with those of several diatoms in the >3 µm size fraction supports the idea that members of this group thrive in the phycosphere and further suggests specific interactions with phytoplankton.
In conclusion, our combined bacterial 16S rRNA and plastid gene metabarcoding approach on size-fractionated samples from DNA and RNA extracts revealed new insights into potential bacteriaphytoplankton associations that could pave the way for further laboratory or in situ studies.In contrast to the 18S rRNA marker gene, the 16S rRNA and rbcL plastid marker genes targeted directly the phytoplankton and gave a complementary picture of the diversity.Furthermore, the more variable rbcL gene was a more informative marker, allowing the discrimination of the dominant phytoplankton species between the two main bloom stations.In analyzing size-fractionated samples, we were able to reveal distinct groups of bacterial taxa  F I G U R E A7 (Continued).
T A B L E A2 Alpha diversity of phytoplankton and bacterial communities in surface depths at the different stations converted into true diversities or "effective species number" (Jost, 2006) from the Shannon Index calculations (DNA).T A B L E A3 Relative abundances of diatoms (primarily order level) at the four stations obtained by microscopy identification (Lasbleiz et al., 2016) of CTD rosette surface samples, compared to those obtained in our >3 µm size fraction samples from 20 m using rbcL marker gene sequencing.T A B L E A5 Proportional abundance of bacterial zOTUs (out of total bacterial zOTUs) correlating to phytoplankton zOTUs and proportional abundance of phytoplankton zOTUs correlating to bacterial zOTUs (out of total phytoplankton zOTUs) for the significantly correlated zOTUs in Figure 6.
sisting of DNA from a synthetic mock community (Mock) of 20 bacterial species containing equimolar (even) rRNA operon counts (HM-782D; Genomic DNA from Microbial Mock Community B, Even, Low Concentration, BEI Resources, Manassas, VA that was amplified as above with the bacterial primer pair.This standard is now obtainable from LGC Standards S.a.r.l.; reference ATCC ® MSA-1002™).Before PCR, DNA samples were diluted in molecular grade water (Merck) to a concentration of 10 ng/µL.We decided to pool the 65-10 µm and 25-3 µm DNA and cDNA size fraction samples for depths 20-150/160 m (see Table 8 µm filter corresponded to a fraction size >0.8 µm <3 µm.DNA/cDNA samples or Mock DNA (1 µL) were amplified in duplicate 10 µL reactions containing 1X KAPA 2 G Fast Ready Mix (Merck) and 0.5 µM of each primer.The PCR cycling conditions were 95°C for 3 min followed by 25 cycles of 95°C for 15 s, 55°C for 15 s and 72°C for 2 s, and a final extension of 72°C for 30 s. Duplicate reactions were pooled and PCR amplification was verified by gel electrophoresis.To normalize the samples before pooling and sequencing, the Sequalprep Normalization Plate (96) kit (Invitrogen) was used according to the manufacturer's instructions.After binding, washing and elution, the environmental and the mock community PCR products were pooled into one tube.The clean-up of the PCR amplicon pool was achieved with the Wizard SV Gel and PCR Clean-Up System (Promega) according to the manufacturer's instructions with elution in 30 µL of molecular biology grade water.DNA was quantified with the Quant-iT PicoGreen dsDNA Assay kit (Invitrogen) according to the manufacturer's instructions.Approximately 700 ng of DNA was pooled with 350 ng of barcoded PCR products from a different project (sequencing the 16S rRNA genes of lichenassociated bacteria) and sent out to the sequencing company Fasteris for library preparation and sequencing (described below).Library preparation involved ligation on PCR using the TruSeq DNA Sample Preparation Kit (Illumina) according to the manufacturer's instructions except that 5 PCR cycles were used instead of 10 cycles.The library was sequenced on one Illumina MiSeq run using the 2 × 300 bp protocol with MiSeq version 3.0 chemistry and the base-calling F I G U R E 1 Location of the four sampling stations with MODIS-Aqua satellite (CLS-CNES) images of surface chlorophyll a (chla) concentrations (µg/L) indicating the presence of the phytoplankton blooms.A3 is the reference bloom station above the Kerguelen Plateau and R is the reference HNLC station.Station F was situated in a bloom above the Polar Front and E5 was located in the Polar Front meander.
Mantel and partial Mantel tests were performed in vegan using mantel() and mantel.partial()based on the Pearson correlation method.Before correlation analysis, environmental variables (depth, temperature, salinity, oxygen, phosphate, silicic acid, nitrite, nitrate, ammonium, and dissolved organic carbon (DOC)) were z-score transformed and a Euclidean distance matrix was calculated.The bacteria and phytoplankton zOTU matrices underwent Hellinger transformation before calculating Bray-Curtis distance matrices.The correlation heatmaps were generated using spls() and cim() in the mixOmics package with clr-transformed data(Chun & Keleş, 2010;Rohart et al., 2017).4| RESULTS AND DISCUSSION4.1 | Environmental contextPhytoplankton and their associated bacterial communities in the particulate fractions (0.8-3 and >3 µm) were characterized at three stations influenced by Fe fertilization to the north-east and southeast the Kerguelen Islands Plateau (stations F, E5, and A3) and at a reference station in the HNLC region west of the plateau (station R)

4. 2 |
Patterns of bacterial and phytoplankton diversity in size-fractionated samples Bacterial and phytoplankton diversity were assessed at the molecular level by generating sequence reads from the 16S rRNA gene (region V3-V4) and the RubisCO large subunit gene (rbcL) respectively from DNA and RNA samples (reverse-transcribed to cDNA).The rbcL gene primers used to target the RubisCO form 1D, present in chromophytic (red) algal groups representing the dominant phytoplankton

4. 3 |
Distinct phytoplankton community structures at the different stationsPhytoplankton communities clustered by a fraction and by station for the upper two surface depths (20-80 m) for both phytoplankton marker genes which showed good agreement apart from the samples from station F (see tanglegram, FigureA2).For both fractions, the communities at R and E5 were most closely related to each other and fell into a cluster with the communities from A3. Whereas the rbcL sequences grouped the 0.8-3 and >3 µm fractions in a separate cluster, the 16S rRNA sequences grouped the A3 and F >3 µm samples with one of the 0.8-3 µm fraction samples.The difference in clustering may be explained by the higher diversity of the rbcL sequences compared to the 16S rRNA sequences which allowed the differentiation of two distinct Thalassiosira sp.zOTUs at A3 and F whereas a single common zOTU was detected in the 16S rRNA data set (Figure2, discussed below).T A B L E 2 Impact of different factors on bacterial community structure tested by PERMANOVA analysis.,E5,R 0.8 µm, >3 µm 20,60-80,150 DNA, RNA DNA(n = 24), RNA (n = 24) 0.24 (0.001) 0.14 (0.001) 0.11 (0.001) 0.09 (0.001) 0.42 Note: R values are indicated with the p Value in brackets.F I G U R E 2 Comparison of phytoplankton taxonomic community structure (order level) at different depths at the four sampling stations in the 0.8-3 µm or >3 µm size fractions as revealed by plastid rbcL and 16S rRNA marker gene sequencing.Diatom orders are shown in shades of green, yellow, and beige from Thalassiosirales to Cymbellales.
. Conversely, in the 0.8-3 µm fraction Prymnesiophyceae largely dominated the R and E5 sequences in the present study.Although Prymnesiophyceae were the most abundant class in the majority of the 0.8-3 µm fraction samples of the bloom stations, attaining >50% of sequences at 80 m at F and 20 m at A3, Bacillariophyceae sequences remained relatively abundant even in this smaller size fraction.At the two surface depths for which rbcL and 16S rRNA plastid sequences were available for comparison (20-80 m, Table diversity and activity of the smaller haptophyte phytoplankton that are difficult to identify by light microscopy.The evaluation of protist diversity at the KEOPS2 stations by 18S rRNA gene sequencing(Georges et al., 2014) revealed the presence of small diatom species at the bloom stations such as Thalassiosira sp., Minidiscus sp. and Coscinodiscus sp. but failed to detect certain species such as Fragilariopsis sp.18S rRNA Haptophyta sequences were also dominated by Phaeocystis antarctica in agreement with our rbcL and 16S rRNA data.However, Bacillariophyceae were underestimated seemingly due to the overrepresentation of alveolates that possess a high rRNA operon copy number(Zhu et al., 2005) and that are known to dominate 18S rRNA gene metabarcode databases(Vaulot et al., 2022).Beta diversity analysis on the rarefied bacterial 16S rRNA gene (DNA and cDNA) sequences revealed that for the bloom stations A3 and F, the community composition was similar between the DNA and RNA (cDNA) samples except for the A3 0.8-3 µm fraction samples at the three depths from 20 to 150 m where the DNA samples clustered separately from the RNA samples.The converse was observed at station R where all DNA and RNA samples clustered separately apart from the 20 m >3 µm sample (FigureA5).Community structure comparisons of marker gene sequences originating from RNA and DNA to reveal activity (protein synthesizing potential;Blazewicz et al., 2013) and characterized by several Proteobacteria families including Alteromonadaceae, Colwelliaceae, Pseudomonadaceae, Halieaceae, Bdellovibrionaceae, Erythrobacteraceae, and Sphingomonadaceae.All the families belonging to the Bacteroidetes phylum were preferentially associated with this larger particle fraction size, as were the Planctomycete families Phycisphaeraceae and Planctomycetaceae and those belonging to the Verrucomicrobia.Low PAN index clades comprised several Oceanospirillales families including SAR86, different SAR11 clades, Moraxellaceae, Marine group B (SAR324), Porticoccaceae (SAR92), and the families Rhodospirillalaceae and Cluster 4 implied potential interactions between Thalassionema, Fragilariopsis, Eucampia, and Naviculaceae with several Flavobacteria zOTUs (Ulvibacter and Fluviicola) and the Gammaproteobacteria Marinicella and OM60[NOR5].Clusters 2 and 3 revealed correlations between the NS9 and NS10 bacteroidetes groups and Pseudoalteromonas with both haptophytes and diatoms (Cluster 2) and between the Gammaproteobacteria OM60[NOR5], the Bacteroidetes Pseudofulvibacter zOTU and the Planctomycete Blastopirellula with diatom zOTUs (Cluster 3).
preferentially co-occurring with the different phytoplankton size classes dominated by Phaeocystis spp.and diatoms, two phytoplankton groups with key ecological and biogeochemical roles.By including metabarcoding sequencing from bacterial 16S rRNA and phytoplankton rbcL transcripts, we were able to add a new dimension to our data by revealing a potential elevated relative activity (RNA:DNA > 1) of all highly correlated Roseobacter zOTUs, Sphingobium, and the Thalassiosira and Chaetoceros zOTUs.The function and importance of the potential interactions we revealed could be investigated in more depth by combining metagenomic and metatranscriptomic analyses of key phytoplankton-bacteria partners in co-culture experiments and natural samples during bloom transitions from diatoms to Phaeocystis.AUTHOR CONTRIBUTIONS Nyree J. West: Conceptualization (supporting); data curation (lead); formal analysis (lead); investigation (equal); methodology (equal); writingoriginal draft (lead); writing-review & editing (equal).Marine Landa: Investigation (equal); methodology (equal), writing-review & editing (equal).Ingrid Obernosterer: Conceptualization (lead); funding acquisition (equal); investigation (equal); methodology (equal); project administration (lead), resources (equal); writing-review & editing (equal).ACKNOWLEDGMENTS The financial support for this research was provided by CNRS-INSU-LEFE-CYBER, ANR-10-BLAN-0614, and the French Polar Institute Paul-Émile Victor (IPEV).Our gratitude extends to Stéphane Blain, the KEOPS2 project coordinator, Bernard Quéguiner, the chief scientist, Captain Bernard Lassiette, and the crew of the R/V Marion Dufresne II for their invaluable assistance and support during the project.We are grateful to the Bio2Mar core facility for granting us access to their equipment.Our thanks go to Karine Leblanc for her contribution to the diatom counts, Yan Liu for her help with data analysis, and Urania Christaki for her insightful discussions.The graphical abstract was created using Biorender.U R E A1 16S rRNA gene (a) and rbcL gene (b) sequencing library sizes obtained from the different sampling sites, depths, and size fractions from DNA and RNA (cDNA).The plastid and bacterial 16S rRNA gene sequences were separated in Phyloseq by subsetting out the Class "Chloroplast," assigned by the Sintax Silva database.F I G U R E A2 Phytoplankton community structure comparisons in the top two surface depths 20 m (1) and 60-80 m (2) at stations A3, F, E5, and R, recovered from two size fractions 0.8-3 µm (0.8) and >3 µm (3) assessed by plastid 16S rRNA and rbcL marker gene sequencing.Dissimilarity dendrograms using Bray-Curtis distances were constructed for each marker gene and their level of congruence was examined by Tanglegram.Sample notation is the Station Name.SizeFraction.Depth.F I G U R E A3 Comparison of phytoplankton taxonomic community structure (class level) at different depths at the four sampling stations in the 0.8-3 µm or >3 µm size fractions as revealed by plastid rbcL and 16S rRNA marker gene sequencing.F I G U R E A4 Comparison of presence and activity of phytoplankton (order level) by rbcL gene sequencing from DNA or RNA (cDNA) extracted from the 0.8-3 µm or >3 µm fractions collected at 20 m at the bloom stations A3 and F. F I G U R E A5 Bacterial community structure comparisons at depths 20 (1), 60-80 (2), 150 (3), and 300 m (4) at stations A3, F, E5, and R, recovered from two size fractions 0.8-3 µm (0.8) and >3 µm (3) assessed by 16S rRNA gene sequencing from DNA and RNA (RNA).Dissimilarity dendrograms were constructed using Bray-Curtis distances.The sample notation is StationName.SizeFraction.Depth for DNA samples and Name.SizeFraction.Depth.RNA for RNA samples.

F
I G U R E A6 Bacterial taxonomic diversity at the order level (DNA) at different depths for the four sampling stations in the 0.8-3 µm (0.8_micron) or >3 µm (3_micron) size fractions.F I G U R E A7 Clustered image heatmaps showing the most significant correlations between the bacterial 16S rRNA zOTUs and phytoplankton plastid 16S rRNA ZOTUs, recovered from DNA (a) or RNA (cDNA) (b).WEST ET AL.|25 of 32 Sampling depths and corresponding environmental parameters are given.The mixed layer depth (Z ML ) is the mean ± sd of all CTD casts performed during the occupation of each station.The Z ML is based on a difference in sigma of 0.02 to the surface value.The euphotic layer is defined as the depth with 1% of surface PAR.
(McMurdie & Holmes, 2013)l.(2015).cDataarefromQuérouéetal. (2015).dValuesareprovidedfor the depth closest to the sampling depth used for sequence analysis within the Z ML (40 m at R and 37 m at A3-2).alignment (v123) and further filtering was carried out to remove zOTUs unassigned at Kingdom level, or those that were not well aligned.The filtered zOTU table, together with the corresponding taxonomy was converted to biom format for import into the R package Phyloseq(McMurdie & Holmes, 2013)for further sequence manipulation and data exploration.
T A B L E A4 Summary of partial mantel tests showing Pearson's R correlations and p Values to explore the abiotic (environmental variables; Env) and biotic (Phytoplankton, Phyto) influence on bacterial community structure in the two size fractions for the top two surface depths (bacterial 16S rRNA and phytoplankton rbcL data set) or all depths (bacterial and plastid 16S rRNA).Tests are presented as the 2 matrices being compared (e.g., Bacteria × Phyto) while controlling for the effects of the third matrix (e.g., |Env).